clc,clear,close;
x=0.99:.0002:1.01;
f1=@(x) x.*x.*x.*x.*x.*x.*x.*x-8.*x.*x.*x.*x.*x.*x.*x+28.*x.*x.*x.*x.*x.*x-56.*x.*x.*x.*x.*x+70.*x.*x.*x.*x-56.*x.*x.*x+28.*x.*x-8.*x+1;
subplot(2,2,1),plot(x,f1(x));
title('f(x)');
f2=@(x) (((((((x-8).*x+28).*x-56).*x+70).*x-56).*x+28).*x-8).*x+1;
subplot(2,2,2),plot(x,f2(x)),axis([0.99 1.01 -2*10^(-14) 4*10^(-14)]);
title('g(x)');
f3=@(x) (x-1).*(x-1).*(x-1).*(x-1).*(x-1).*(x-1).*(x-1).*(x-1);
subplot(2,2,3),plot(x,f3(x)),axis([0.99 1.01 -2*10^(-14) 4*10^(-14)]);
title('h(x)');
subplot(2,2,4),plot(x,f1(x));
hold on;
plot(x,f2(x));
hold on;
plot(x,f3(x));
title('三者对比');
hold off;